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Abstract 

An analogue of Kolmogorov's superconvergent perturbation theory in classical 
mechanics is constructed for self adjoint operators. It is different from the usual 
Rayleigh-Schrodinger perturbation theory and yields expansions for eigenvalues and 
eigenvectors in terms of functions of the perturbation parameter. 

PACS Code: 03.65.-W, 31.15+q, 02.30.Mv, 02.90. +p 

1. Introduction. About a century ago Lindstedt Poincare 0, and von Zeipel [0 
developed a perturbation theory for Hamiltonian systems in classical mechanics (CM). 
The method was mainly used in celestial mechanics but in most cases failed to converge 
due to the appearance of small divisors. It was only in the fifties that Kolmogorov 
proposed a new quadratically convergent perturbation method for Hamiltonian systems 
in CM and this (nowadays called superconvergent) method was essential in Arnold's proof 
of the KAM theorem [§. 

In this letter an exact analogue of Kolmogorov's superconvergent method is con- 
structed for self adjoint operators. So far "superconvergent" is only a name for this new 
method and a detailed functional analytic investigation will have to determine how much 
the new perturbation theory constructed here is an improvement on existing schemes. 
There are, however, indications that this will be so: The first comes from the fact that we 
have also constructed an analogue of the Poincare-von Zeipel perturbation theory for self 
adjoint operators and have shown that this analogue is identical to the usual Rayleigh- 
Schrodinger perturbation theory 0. Since in CM the superconvergent method is a vast 
improvement on the Poincare-von Zeipel method we may expect the same in Quantum 
Mechanics (QM). The second comes from initial numerical studies. 
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2. General Algorithm of the Superconvergent Method. Let Hq be the unperturbed 
Hamiltonian (CM: function on phase space, QM: operator) and let 

OO p 

p=o 

be the perturbed Hamiltonian where the do not depend on the perturbation parameter 
e. Moreover, let 

OO p 

W-{e) := E -W;^, (2) 
p=o t'- 

and let —W{e) be generators (CM: functions on phase space, QM: operators) of (CM: 
canonical, QM: unitary) flows {<^^{e))^^ with "time" e. This means that the transforma- 
tions $'^(e) satisfy the initial value problem 

= o$"(e)*, $"(0) = id (3) 

ae 

where in CM the $"(e)* act on any phase space function A, via 

$"(e)*(A) := Ao$"(e) (4) 

and where in CM 

adW-{e){A) A}, (5) 

and {■, ■} denotes the Poisson bracket. In QM self adjoint operators —W^{e) generate 
one parameter groups of unitary transformations ($"(e))~^ such that the $"(e) solve the 
same initial value problem as given in (|^) but where the $"(e)* act on operators A via 

$"(e)*(A) := <l>"(e)"i A $"(e) (6) 

and where now 

adW''{e){A):=^[W'^{e),A] (7) 

and [■, ■] denotes the commutator. With the understanding that the appropriate definition 
(^ or (depending on whether one is dealing with CM or QM) is chosen we proceed to 
present Kolmogorov's superconvergent method. In the followng it will be useful to have 
an expansion for $"■(£)* in terms of operators independent of e. Writing 

OO p 

= E -,t; (8) 

p=o P- 

one finds that the are recursively determined by Tq = id and 



(=0 
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In CM the are differential operators such that Tp{A) means apphcation of to the 
phase space function A whereas in QM they are operators acting on operators. 

Kolmogorov's method consists of finding the generators VF"'(e) in such a way that 



K"~i(e) := $"-i(e)* o $"-2(e)* o ■ ■ • o $i(e)*(i/°(e)) 



(10) 



has an expansion in e 



oo 



p=0 t'- 



fill 



with "integrable" terms up to order 2"^^ — 1 (whereas in CM the term integrable has a 
well defined meaning this fails to be the case in QM; however, we only use the notion of 
integrability to motivate the super convergent method, all its equations are well defined in 
CM without any need to define integrability and will be well defined operator equations 
in QM even though the term integrability is not). Then a new perturbed Hamiltonian 



(e) := K"^ [e) is defined such that its unperturbed part i/g" consists of the in- 
tegrable part of X""^ and its perturbation is of the order 2"^^. We can summarize the 
iterative procedure as follows (for a more detailed discussion we refer the reader to [0 
with the warning that our notation differs from theirs): 



result of 
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choice of 



result of 



(int. up to 0(2 



n-l 



1)) 



P 



0, 



1 < p < 2 



71-1 



adW;{H^ 



n— 1 > 



H'^-^ = 0, 1 < p < 2"-i 
(no perturbation) 



iY^-i - 2"-i < p < 2" 



with adH^-^i 



wi 



0, 



2" < p 



leads to 



P 
P 



0^ 

zjn—l 

Hp 



l<p< 2"-i 

2n-l < p < 2" 



= ^0 

(int. up to 0(2" - 1)) 



= 0, 1 < p < 2^ 
(no perturbation) 



TTn 

Hp 



2" < p 



2" < p 

The meaning of ~ is to be understood as follows. Assume H^~^ is of the form given in the 
left column and choose Wp = for 1 < p < 2"~^. Then one finds first K^ 
0, for 1 < p < 2"-i and 



'■0 



Zjn-1 TT-n 



P 



adW;{Hr') + H'p 



n-l 



2""^ < p < 2". 



:i2) 



Since Hq is already integrable up to 0(2" — 1) and $" should improve this, one would 
like to have that the Kp for 2"~^ < p < 2" are integrable. Consequently the crucial point 
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in this procedure becomes the construction of the and some ^ for 2" ^ < p < 2^ 
such that 

adW;{H^-^) = - H;-^ (13) 



where H^~^ is such that 



adH;-\m-^) = Q. (14) 



Then it follows from (|T2]) that K'^ = H^~^ and from (|Tj) that they commute with the 
unperturbed Hamiltonian Hq~^ of the previous step. The table summarizing the method 
omits these intermediate steps and shows only the results once (|13D and (|l^ have been 
solved. 

In CM (|13|) and (p!4D can be satisfied with the help of the averaging method and 



Hp ^ turns out to be the average of over the angle variables of the tori which are 

assumed to be compact. It should be noted, however, that the Hp~^ in general depend 
on e such that (|T3p is actually an expansion in terms of functions of the perturbation 
parameter. Since Hq~^ is assumed to be integrable (|14D implies that Hp~^ and hence Hq 
are integrable as well. 

Hence, after the nth transformation one has a Hamiltonian H"'(e) where 




=2^-1 P 



is integrable and the perturbations are of order 2" or higher. 

3. Superconvergent Method for Self Adjoint Operators. With our previous notation 
equations (|T3p and ( [T^ ) become operator equations in QM which have to be solved for each 
step. We can construct operators Wp, Hp~^ satisfying (|13D and ([I^) with the help of the 



following quantum analogue of the classical averaging procedure P] which is a modified 
version of an idea used by Weinstein in the context of pseudodifferential operators: 
Let A, B be self adjoint operators such that 

B^^\t) := exp(-^tA) 5 exp(^tA) (16) 
n h 

B^^^ := \im^ C B^^\t)dt (17) 







S^^\B) := lim - r dt t ds (b^^\s)-B^^A 

T^oo T Jo Jo \ J 



exist and 



lim = (19) 



then the following holds 



adB'-^^A) = '-[B^^\a]=Q and (20) 

ad S^^^ (A) = Us^'^MB), A] = B^^^ - B. (21) 
n 
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Proof: 



adB^^'^ (A) = lim 1 UB^^\t) , A]dt = hm i ^B^^\t)dt 

T^oo 1 Jo h T~*oo 1 Jo at 



= lim = 

T^oo T 

by assumption and thus 

T^oo 1 Jo ^ ' 

which completes the proof. Setting now = S^^o solves (jlSj) and (0) in 

QM with Hp^^ = Hp~^ ° =: Hp where the last equation introduces a simplified 
notation. 

4-. Interpretation of the Algorithm for Self Adjoint Operators. The results of section 3 
show that we can execute Kolmogorov's super convergent perturbation algorithm in QM 
as well. But what has been gained? By construction H"'{e) and the original perturbed 
Hamiltonian H^{e) are unitarily equivalent: 

i/"(e) = $"(e)-i ■ • • <^\e)-^ H\e) $^(e) • ■ ■ <l>"(e) (22) 

and thus have identical spectra. Since, after the first transformation [Hq,Hi^^^] = 0, 
we can diagonalize Hq and Hi^^^ simultanously which permits us to diagonalize Hq := 
Hq + eHi^\ After the second transformation we have 

[HIh}'^]=Q = [HIh}'^] (23) 

and thus we can diagonalize 

ff„= = i/S + ei?7'»' + f7?;'" + ^7^'" (24) 

etc. In this way we arrive after n transformations at a diagonal i^o (e) whose eigenvalues 
E^[e) approximate those of H"'{e) and thus of H^{e). The eigenvectors of are 
known by construction and they approximate those of if". Because of ( P^ it follows 
that <l>^(e) ■ ■ ■ $"(e)|j)"' are the appropriately approximated eigenvectors of the original 
perturbed H^{e). In formulae: If H^{e) \j){e) = Ej{e) \j){e) then 

E,{e) = E]ie)+0{r) (25) 
|j)(e) = $i(e)-..$"(e)|jr + 0(2'^). (26) 
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Since the Hp depend on e the E^{e) are no longer a pure power series expansion but 
should be viewed as an expansion in terms of functions of e. These functions themselves 
have in general infinite scries expansions in e so that already low orders i?"(e) contain 
contributions of all orders in e. The example treated in the next section will make this 
point clearer. 

5. An Example: Discrete and Nondegenerate Spectrum. It remains to show that the 
method constructed her is truly distinct from the usual Rayleigh-Schrodinger perturbation 
theory. For this purpose we sketch the results for the case of a Hamiltonian Hq = 
J2j \j) Ej {j\ with purely discrete and nondegenerate spectrum and a perturbation only 
linear in e, i. e. H^ = J2j,k \j) ^jk {k\ and H^ = for p > 2. In this case the eigenvectors 

of Hq = Hq + eHi^^ and \j) of H^ coincide because Hq and Hi'^^ commute and Hq 
has nondegenerate spectrum. Continuing this argument inductively it is easy to see that 
Hq and all Hq have the same eigenvectors \j). By straightforward calculation one finds 

- 7Eb)|^^(^l (27) 

= ^? + E E nEf\ (29) 

where Ef' := <l — 1 and the E^ are the eigenvalues of H'^. From 

the algorithm presented in section 2 one obtains 

Hi = \WlHi'^ + ii'^ (30) 
Hi = Q\w^,[WI,H;^'^ + 2H',]] (31) 

Hi = Q\wl[Wl[W^,Th^'^ + 3H',]]] + 3^[WlH-,^'^+Hl]. (32) 

The appearance of the denominators E^~^ — E'^~^ in leads to a perturbation the- 
ory substantially different from the usual Rayleigh-Schrodinger theory. However, this 
difference only shows up in the fourth and higher order terms since in Hi and Hi only 
commutators with Wl appear which contains only the usual Rayleigh-Schrodinger de- 
nominator Ej — E^. Indeed, one finds 

E^'^ = V,, (33) 
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1 7-i(3) VjkMkmVmj \^jk\ ^jj fQK\ 

showing that up to 0(3) the new method coincides with the usual Rayleigh-Schrodinger 
perturbation theory. But iJ| contains W2 with denominators of the form Ej — E\ = 
Ej — E^ + eiVjj — Vkk) which are functions of e and which appear as denominators in 



{E^ - E^Ej - El) 

VjlVlkVkmVmj /I 1 \ / 1 



+ 6 E 



j^l^k^j^m^k ^k ^k) \^k ^ra , 



+ 12 E Vji^ikVk 



Vu-V, 



kk 



kj 

j=a^k^j L v-^i ^u\-^i ^kj\-^k 

Vu - V,, I 1 



{E]-E^){E^,-Et){El-Ef) 
\ 



[E^ - Ef){E] - El) \Ef - El El - El 

\k ~ Vjj 



/ y-^l ^k ^k ^] 

^ Vkk-V,, ( 1 1 



{E] - El){E] - El) \E^ - Ef El - El J 

^ wS.^/^'^'^'™^'"'' I (^^ - E^iE^ - El) ~ W^k, 

3 (I 1 M 

+ 



{El - El){El - El) \E^ - El El - El 



As an example we quote the result for the ground state correction up to £'0'^'' for the 

cP 



quartic perturbation = of the harmonic oscillator Hq = — + a;^ (with ^ = 1) 



0^ ^4 16 64 



3(1317760 + 12935472e + 36433368e2 + 25183305e 
2048(4 + 9e)(4 + 15e)(4 + 21e) 



3\ 



-e' + 0(6^) 



which approximates the numerically computed eigenvalues much better than the stan- 
dard fourth order Rayleigh-Schrodinger correction 

E^Je)ns = 1 + -e - -e^ + — - ^^^^e^ + 0{e'). (36) 
ov )Rb T- ^ 16 64 1024 ^ ^ > ^ > 

This example also illustrates the statements made earlier that the lower order corrections 
E^l'^ already contain infinite power series of e which may be the reason for improved 
convergence. 

6. Conclusion. It should be emphasized that the method presented here is very 
general and in principle applicable to any self adjoint operator A and perturbation B 
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provided they (and the higher order operators) satisfy the conditions needed to construct 

B^^^ and S^^^{B) in ( ^n\ ) and ([TB|). It is also designed to include analytic perturbations 
which are not necessarily only linear in the perturbation parameter. Moreover, it gives 
the corrections as integrals of the form B which, if desired, can be exhibited as sums 
over intermediate states as was done in section 5, but which may be evaluated directly 
circumventing the calculation of these sums which may sometimes be impossible. This was 
already shown to be an advantage of the formulation of the usual Rayleigh-Schrodinger 
theory as an analogue of the classical Poincare-von Zeipel perturbation method in . 



The method will be presented in more detail in a forthcoming longer paper |10 
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